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Abstract 

A thermo-viscoplastic computational meth- 
od for hypersonic structures is presented. The 
method employs a unified viscoplastic constitu- 
tive model implemented in a finite element ap- 
proach for quasi-static thermal-structural anal- 
ysis. Applications of the approach to convcc- 
tivcly cooled hypersonic structures illustrate the 
effectiveness of the approach and provide insight 
into the transient inelastic structural behavior 
at elevated temperatures. 

INTRODUCTION 

The commitment to develop the National Aerospace 
Plane has generated resurgent interest in the technol- 
ogy required to design structures for hypersonic flight. 
Such structures will be exposed to aerodynamic heating 
of unprecedented magnitudes. As the vehicle accelerates 
or decelerates at hypersonic speeds in the atmosphere, 
shocks will sweep across the vehicle and interact with 
local shocks and boundary layers. These interactions in- 
troduce severe local pressures and heating rales. A recent 
experimental study (ref. 1) of interacting shock waves on 
a cylindrical leading edge shows heating rates ten times 
undisturbed levels. 

Leading edges of engine structures (Fig. 1) present a 
significant design problem because of intense local heat- 
ing and pressures. Analysis of the fiow, thermal and 
structural behavior present serious computational chal- 
lenges to analysts because of the inherent nonlinearilies 
in all aspects of the multidisciplinary problems. Some 
of the critical computational issues are identified in ref- 
erence 2. Critical issues include the difficulties involved 
in (1) analyzing the viscous, compressible (low and pre- 
dicting the high local aerodynamic heating, (2) model- 
ing and analyzing multimode unsteady heat transfer in a 
high temperature convcctivcly-coolcd structure, and (3) 
simulating the transient, nonlinear thermal-structural re- 
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sponse for rapid temperature changes. Preliminary struc- 
tural analysis of an impingement cooled leading edge (ref. 
3) showed high local plasticity that seriously degraded 
the structure’s load carrying capacity at elevated tem- 
peratures. A recent thermostructural analysis with ex- 
perimental verification (ref. 4) of cowl lip designs con- 
firmed that inelastic effects occur and can be significant. 
In the experimental study, two specimens failed due to 
burn-through because of intense local heating or because 
of loss of cooling. 

The purpose of this paper is to present a thermo -vis- 
coplastic computaliontal method for hypersonic struc- 
tures subjected to severe local unsteady healing. The 
analysis employs a unified viscoplaslic constitutive model 
implemented in a finite element approach capable of pre- 
dicting rate-dependent plasticity effects for temperatures 
up to about 75 percent of the melting point. Rate- 
dependent plasticity effects are known to be important 
at elevated temperatures. 

Unified viscoplastic constitutive models have evolved 
over the last twenty years to provide a means for ana- 
lytically representing a materials response from the elas- 
tic through the plastic range including strain rate de- 
pendent plastic fiow, creep and stress relaxation. The 
theories are guided by physical considerations including 
dislocation dynamics and are based on the principles of 
continuum mechanics. The first multidimensional formu- 
lations of elaslic-viscoplastic constitutive equations was 
due to Bodner and Parlom. Since then a number of 
constitutive models have appeared; many of these the- 
ories are summarized in review articles that appear in 
reference 5. A NASA -Lewis sponsored research program 
(HOST) conducted by the Southwest Research Institute 
recently concluded a four year research effort (ref. 0 -7) to 
further develop unified constitutive models for isotropic 
materials and to demonstrate their usefulness for analysis 
of high temperature gas turbine engines. One result of 
this study is material property data for high temperature 
nickel-based alloys over a wide temperature range. The 
unifed models employed were those of Bodner-Partom 
and Walker. 

Unified viscoplastic theories have been implemented 
by a number of finite element researchers. Detailed stud- 



ies of several rate-dependent plasticity models and tlieir 
numerical implementation using adaptive finite element 
methods were conducted by Bass [8j and Bass and Oden 
(9, 10). Under the NASA HOST program, the Walker 
model was implemented in the MARC finite element 
program (ref. 11) and used to analyze the thermo- 
viscoplastic response of a turbine blade under simulated 
dight conditions. In another recent finite element appli- 
cation (ref. 12), the Bodner-Parlom and Walker theories 
were compared for a thin circular plate subject to highly 
localized, transient heating. 

In this paper the Bodner-Partom constitutive model 
is employed, and the finite element approach developed 
in references 8-10 for the isothermal case is extended to 
include thermal effects. The paper begins by describ- 
ing the finite element approach used for thermal anal- 
ysis of convectively cooled hypersonic structures. Then 
the thermo-viscoplastic structural analysis is presented. 
The initial value problem is formulated, the constitutive 
model is described and the finite element viscoplastic so- 
lution method is presented. Next, applications of the 
approach arc made to convectively cooled structures. An 
appendix presents thermal and structural data for the 
convectively cooled hypersonic structure used in the nu- 
merical computations. 

THERMAL ANALYSIS OF HYPERSONIC 
STRUCTURES 

Convectively cooled structures are strong candidates 
for use in hypersonic flight vehicles. For hypersonic flight, 
some leading edges and panels require active cooling sys- 
tems to keep structural temperatures within acceptable 
ranges. The internal (low in the coolant passage has a 
predominant role in the thermal response of a hypersonic 
structure subject to external heating. A cross-section of 
a typical convectively cooled structure is shown in Fig. 2. 
An aerodynamic skin and a coolant passage with internal 
heat exchanger protect the primary structure from the 
aerodynamic heating. The thin, typically metallic, aero- 
dynamic skin transfers the energy of the aerodynamic 
heating to a low temperature coolant flow through the 
heat exchanger fins that connect the aerodynamic skin 
to the primary structure. In a typical engine structure, 
the coolant is cold hydrogen that later is used as the 
propulsion system fuel. 

Heat transfer in the aerodynamic skin consists of con- 
duction combined with surface radiation. Heat transfer 
between the aerodynamic skin, the heat exchanger sur- 
faces and the primary structure is by conduction at the 
solid -fluid interface. The finite clement representation of 
conduction heat transfer with radiation boundary con- 
ditions follows the standard procedures described in ref. 
13. 

The dominant mode of heat transfer in the coolant 
(low is forced convection. The representation of the heat 
transfer in the coolant passage is the critical step in the 
heat transfer analysis. There are two basic representa- 
tions that can be employed. The first, denoted here as 
the engineering model, is based upon a number of as- 
sumptions that greatly simplify the problem into a single 
energy equation with a specified mass flow rate. De- 
tailed computation of the fluid velocity components and 


temperatures is not required. The second representa- 
tion idealizes the coolant (low as a continuum model, and 
the partial differential equations describing conservation 
of mass, momentum and energy are solved simultane- 
ously to obtain fluid velocity and tmperature distribu- 
tions. This latter model is the most accurate, but it is 
also considerably more expensive than the engineering 
model. In this paper the engineering heat transfer model 
is employed. 


Engineering Model for Coolant Passages 

The basic features of the engineering model are devel- 
oped from the idealization shown in Fig. 3. A segment of 
the coolant passage of width to is shown; for simplicity, 
only the upper one-half of the coolant passage with the 
aerodynamic skin is shown. The engineering formulation 
(ref. M) is based on the following assumptions: 

1. The thermal energy state of the fluid is character- 
ized by the fluid bulk temperature Tp which varies 
only in the flow direction, i.e., Tp(x,t). 

2. The (low is represented by the mass (low rate lii 
in the coolant passage specified by m = ftp ApVp 
where pp is the coolant density, Ap is the cross- 
sectional area of the coolant passage, and Vp is the 
coolant mean flow velocity. 

3. A convection coefficient /i is defined such that the 
heat flux qp transferred between the structure and 
the coolant may be expressed as 

qp = h(T s - Tp) 

where Tg(x,t) denotes the structural temperature 
at the fluid-solid interface. 

4. The convection coefficient h may be expressed as 
a function of the fluid bulk temperature alone by 
using analytical/empirical equations for the Nusselt 
number, 


where D is the hydraulic diameter of the coolant 
passage, and kp is the thermal conductivity of the 
fluid coolant. The convection coefficient as well as 
the solid and fluid thermal parameters may, in gen- 
eral, be temperature dependent. 

With these assumptions, energy balances on the aerody- 
namic skin and coolant give the governing conservation 
equations. 

Fluid: 

[ krAr ^t] - u M T S-Tp)+PpCpAp^ =0(1) 


Solid: 

ksA s %jf) + wh(T s - Tp) + acV 

+ PscsAs^jjf- = q 


on 



where the subscripts F and S denote the fluid and solid, 
respectively. In these equations, c denotes specific heat, 
a is the Stefan- Boltzmann constant, and c is the surface 
emissivity. Because of the temperature dependence of 
the thermal parameters and the radiation term, Eqs. (1) 
and (2) constitute a nonlinear set of partial differential 
equations. 

The heat exchanger fins are not included explicitly in 
this model. However, the heat transfer between the heat 
exchanger fins and coolant can approximately be taken 
into consideration through the use of an effective width 
w which represents the area over which the convective 
heat exchange occurs. 

Finite Element Formulation for Coolant Passages 

A typical finite element representing Eqs. (1) and (2) 
is characterized by fluid and fluid-solid interface nodes. 
The element shown in Fig. 4 has two fluid nodes (I 
and J), and two fluid -solid interface nodes (ft and L). 
Within the clement, the fluid and solid temperatures are 
expressed as 


7> = [/v(x)|{r,.} 


Ts = [N(x)\[T s ) (3) 


where [/V(x)] are the element interpolation functions. 
Following usual finite element procedures (ref. 13), the 
discretized equations for an element of length L may be 
derived in the form 
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where the clement capacitance matrices arc given by 


l c,,\ = j L p,,c F A F (N}[N}dx 

JO 

[Cs] = j\sc s As{N)[N\dx 


(5) 


and the element conductance matrices are 


The element equations given in Eqs. (4) show that the 
coolant passage model can be regarded as an assembly 
of elements where each element represents a single heat 
transfer mode. Thus we can represent the coolant pas- 
sage in terms of two convective elements: (1) a mass- 
transport element (Fig. 5a), and (2) a surface convection 
element (Fig. 5b). The mass transport element repre- 
sents the downstream convective heat transfer due to 
mass flow; it has a capacitance matrix [CV-] as well as 
conductance matrices [A'„] and [ftp], The surface con- 
vection element represents the convection heat exchange 
between the coolant fluid and the solid; it does not have a 
capacitance matrix, but it has a conductance matrix [A'*) 
linking the solid and fluid nodes. These coolant passage 
elements are assembled into the finite element thermal 
model for the complete convectively cooled structure that 
has conduction elements, radiation elements, and surface 
elements for convection to a specified convective exchange 
temperature. 

The unsteady thermal analysis is nonlinear because of 
temperature dependent thermal properties and surface 
radiation. The equations are solved by time marching 
with the Crank-Nicolson algorithm; at each lime step, 
the nonlinear algebraic equations are solved by Newton- 
Raphson iteration. 

TIIERMO-VISCOPLASTIC STRUCTURAL 
ANALYSIS 

The behavior of a thermo-viscoplastic structure sub- 
jected to aerodynamic heating is analyzed assuming that: 
(l) thermo- mechanical coupling in the conservation of 
energy equation can be neglected, (2) the structural re- 
sponse is quasi-static, and (3) deformations are infinites- 
imal. With these assumptions, an unsteady thermal 
analysis may be performed first to determine the tem- 
peratures. Then, using these temperatures, the struc- 
ture’s viscoplastic response is determined. The solution 
is thus obtained by separately solving initial boundary 
value problems for first the thermal and then the struc- 
tural response. 

Initial Value Viscoplasticity Problem 

Consider a viscoplastic structure occupying a region 
fl with boundary dtt. The behavior of the structure is de- 
scribed by the following system of differential equations: 

1. Equilibrium in rate form, 


r L dN 

[K«] = J a «»cjr{Af}[— \dx 

[ft),] = f iuh{N }[ft]dx 

Jo 


°<j.i + i>i = 0 (7) 

where <r,y denote components of the stress tensor, 
6, are the body force components per unit volume, 
and the summation convention is employed. 
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( 6 ) 


2. Kinematic relation for velocity gradients, 

+i Z = + M W 

where t,y denotes the total strain components and 
superscripts E and P denote elastic and inelastic 
strain components, respectively. The components 
of the displacement rates arc ti,. 


3. Constitutive relations 

g I j = EijkiijH — Eiju<*uAT 

£,'■ = fiA a i>i Z k) (no sum) (9) 

Z; = 9<(<?>,, z k) (no sum) 

wliere E,,u re|)resenls Hooke’s tensor of elasticity 
parameters, r»n arc coinoncnts of a tensor of ther- 
mal expansion parameters, and AT represents the 
rate of the change in temperature. Doth and 
on are temperature dependent. The constitutive 
functions are /q and gt wliere Z* represents inter- 
nal state variables. These functions and state vari- 
ables characterize the viscoplastic response of the 
material. 


The description of the problem is completed by pre- 
scribing the boundary and initial conditions, 


it = u, on d(l\ 

ffijiij = a, on an, 


(10) 


where tl, are prescribed surface displacement rates, n / 
are the components of a unit normal vector, and Oi are 
prescribed surface traction rates. The initial conditions 
include specifying the displacements, stresses and inter- 
nal state variables, i.e., 


u,( *,o),(T 0 (*,o),z,(*,o) * € n. 


Constitutive Model 


The Uodner-Parloin constitutive model is of the inter- 
nal stale variable type that is based on phenomenological 
observations and supported by physical concepts related 
to dislocation dynamics. The model has gone through 
several modifications and was extended for anisotropic 
work hardening materials. The current model (ref. 5-6) 
also includes temperature effects. 

1 . Flow Law, 

For the inelastic strain rate component, the isotro- 
pic form of the Prandtl- Retiss law is assumed 


= AS,, 


4 = 0 A > 0 


(H) 


where S,j are the deviatoric stress components 
given by 

Sij = O ij — ~8ij(Tkk (12) 

and = 0 denotes plastic incompressibility. 


2. Kinematic Equations, 
Squaring Eq.(U) leads to 


A J = Df/J, (13) 

where l)£ and J, are the second invariants of the 


strain rale and stress tensors, respectively. For ex- 
ample, the deviatoric stress invariant is 

Ji — -S^Stj (14) 

The relation governing inelastic deformations is the 
“kinetic equation”, and the form taken by Bodncr- 
Partom is 


= D\ exp[-(ZV3^) n ) (15) 


where D 0 is the limiting strain rate in shear, n is a 
temperature-dependent material parameter, and Z 
is interpreted as a load history dependent param- 
eter (herein called the internal state variable) that 
represents the hardness of the material with respect 
to resistance to plastic flow. Combining Eqs.(li), 
(12), and (13) gives 



3. Evolution Equations of Internal State Variable, 

The internal state variable Z consists of isotropic, 
and directional components, 

Z = Z' + Z D (17) 


The evolution equation proposed for the isotropic 
hardening component (ref. 5) is 


Z'(t) = 


m,|Z, -Z,(0]H , „(Z)->l,Z, 


z'm - z* 


r I 


(18) 


with the initial condition, Z'( 0) = Z 0 - In the first 
term, Z, is the limiting (saturation) value of Z 1 , 
mi is the hardening rate, and the plastic work rate 

IKp = (19) 

which is taken as the measure of hardening. Z, is 
the minimum value of Z f at a given temperature, 
and /li and r, arc temperature dependent material 
constants. 

The evolution form of the directional hardening 
component (ref. 5) is defined as 

z D (i) = /M<K(<) (20) 

where u tJ are the direction cosines of the current 

stress state, 

xi, j = <r„(f)/(or w (f)<r w (t))^ (21) 


The evolution equation for has the same gen- 

eral form as that for isotropic hardening but has 
tensorial character, 


/mo = »*,[z 3 Mii(o -ftj(o]iiyo 
I l/MO/MCll 


— AiZ, 


1 


MO 


( 22 ) 


where 

•fi(0 = /M0/IA«(0A|(0I* (2.1) 

and 

/MO) = 0 (2d) 

As in Eq.(4.12), m 3 is the hardening rate. Ai and 
r 2 are temperature dependent material constants. 

Other variants of the isotropic and directional hardening 
variables arc described in refs. (5-7). These references 
also present data for the parameters used in constitutive 
models for several materials. The particular values used 
in the present study are presented in the Appendix. 

Finite Element Formulation 

The finite element formulation follows the approach 
of Ref. 9, and the details of the weak formulation are 
presented there. Since the approach presented in Ref. 9 
is for the isothermal case, this formulation extends the 
previous approach by the inclusion of temperature effects. 

The finite element approach approximates displace- 
ment rates within an element by taking the displacement 
rates as 

{d} = [/V](6} (25) 

where [/V] are the interpolation functions, and {5} rep- 
resent the nodal displacement rates. Using the strain- 
displacement equations in rate form, Eq. 8, an element’s 
strain rates can be computed as 

(r) = [B]{i) ( 26 ) 


where [/}] is the strain-displacement matrix. Following 
usual finite element procedures, the finite clement equa- 
tions for a typical element are form as, 


[A'(T)](6} = 

- {F,) + {Ft} + {K) + {Fb) 

(27) 

where [A’(7’)] is the clement stifTness matrix, and the 
terms on the right-hand side are element load vectors 
due to the rate of plastic strains, temperature, surface 
tractions and body forces, respectively. These matrices 
are defined by 

|A'(T)) 

= J n \R\ T [E(T)\[B\dn 

(28) 

{?,) 

= [ (D| T (E('r)|{€' , }dn 
J ft* 

(29) 

(A) 

= / (D) T [F('r))(«(T))Afdn 

./f). 

(30) 

{^} 

= f [N] T {if}Js 
Jdn e 

(31) 

{^b} 

= I [/v] T (i>}dfi 

yn« 

(32) 


where fi e denotes the clement volume, and dU c denotes 
an clement surface where tractions are defined. 

The temperature affects the viscoplastic structural 
analysis directly in three ways: (I) the elasticity matrix 
{ Z?( T’)) and the coefficients of thermal expansion (o(7')) 


depend on temperature, (2) nodal loads (Fj) depend on 
the local temperature rales, and (3) several parameters in 
the Bodner-Partom constitutive model are temperature 
dependent. 

Viscoplastic Solution Method 

Since the present approach uses an uncoupled (some- 
times called one-way coupled) formulation, the thermal 
problem is solved first followed by the viscoplastic anal- 
ysis. The transient thermal problem is solved by time 
marching with a time step A It, and the nodal temper- 
atures at successive times are obtained. These 

temperature vectors are used as input to the structural 
analysis. 

The first computation in the structural analysis is to 
solve an initial statics problem if the initial temperature 
distribution T(x,0) is not equal to a uniform reference 
temperature or if any initial static loads are present. The 
results of this analysis are the initial conditions (displace- 
ments and slresses) for the transient viscoplastic analysis. 

The viscoplastic analysis time-marches with a time 
step At,. Experience has shown the time step required 
for the structural analysis is usually smaller than for the 
thermal analysis, i.e. At, < Afj-. At intermediate times 
in the structural analysis, the temperatures arc linearly 
interpolated from the temperatures known at the begin- 
ning and end of the larger thermal time intervals. 

The strategy employed in the viscoplastic algorithm 
is as follows: with the initial distribution of stress, tem- 
perature and internal variables specified use the equilib- 
rium condition (Eq. 27) to obtain the nodal displacement 
rates. Then integrate the constitutive equations forward 
in time at the element Gauss integration points. With 
updated values of the stress, temperature and internal 
variables at the new time, the equilibrium equation is 
solved again. This sequence of determining the nodal 
displacement rates, then advancing the constitutive equa- 
tions in time is continued until the desired history of the 
initial boundary-value problem has been obtained. 

Thus, the algorithm proceeds through the following 
steps: 

1. At time t, initialize <7,y, Z, for each element; 

2. Calculate e£ = /,,(cr;y , Zj. ) for each element; 

3. Assemble and solve (A )(6) = (F) ; 

4. Calculate e,y for each element, { c) = (D)(6); 

5. Calculate d , t for each element, (d) = (F) (t — c P } — 
[F)aA7' ; 

6. Calculate Zy for each element, Z, = jy(<T,y,Z*); 

7. Integrate dyy, Z forward for each clement to get <j tJ 
and Z, at t + A/,; 

8. If t + Af, < go to 2, otherwise stop. 

The computational method above has been presented 
for a constant time step At,. Computational experience 
by several investigators (see Ref. 7-12) indicates that a 
very small time step can be required because of the “stiff’ 
nature of the ordinary differential equations describing 


the internal state variables. To gain improved cflicioney 
and reliability a variable time step algorithm has been 
implemented. 

Variable Time Step Algorithm 

In the constant time step viscoplastic algorithm the 
internal stale variables are advanced in time with the con- 
ditionally stable Enler forward dilference algorithm. The 
variable time step algorithm is a modified Euler scheme 
using a truncation error criterion (ref. 15) to adjust the 
time step. 

For simplicity, consider the single ordinary differential 
equation, 

y = f{yj) ( 33 ) 

The solution is advanced using a predictor -corrector 
scheme. The predictor phase consists of an Euler step: 


1/f+ai = >Jt + 

(34) 

rfU, = + ^0 

(35) 


An error indicator E ( Ref. 15) is then computed from 

E = |A<( ^V ~ -- 1 (36) 

^\yi+At I 

The error indicator is next compared with a preset error 
criterion, and if the criterion is met, the time step is suf- 
ficiently small enough to proceed to the corrector stage. 
Otherwise, the predictor phase Eqs. (34-35) is repeated 
with a smaller time step. 

The corrector phase is the modified Newton scheme, 
y*v, = (lit + !/f + ai )/ 2 ( 37 ) 

i/?+<u = 'Jt + kly ava (38) 

A flow chart depicting the adaptive scheme is shown 
in Fig. 6. The flow chart shows how the time step is 
either reduced or increased depending on the error indi- 
cator, Eq. (36). The flow chart shows that the time step 
is reduced or increased by a factor of two. This approach 
is effective, but an alternate scheme also has been used 
where the new time step is based on the error E. In this 
scheme 

A l = A t ] j^ (39) 

where E opt = \/E m „ % E mm . Error tolerances and 

E mm are specified by the user, typical values being 
= 0.05, E„„„ = 0.005. 

APLICATIONS 

As mentioned in the introduction, scramjet engine 
structures can experience intense local healing due to 
shock interactions. Two models of such structures are an- 
alyzed for heat loads representative of shock interactions. 
Material properties used (see the appendix) represent a 
high-temperature nickel-based super alloy. Computa- 


tions are performed first for a simplified one dimensional 
model to provide insight into the computational proce- 
dure. Then a two-dimensional model of a realistic con- 
vectively cooled structure is analyzed in detail. 

Simplified One-Dimensional Model 

To gain a preliminary understanding of the behavior 
of a convectively cooled structure subjected to convec- 
tive heating, a simple ID bar model was investigated. A 
segment (Fig. 7a) of the aerodynamic skin is modeled 
assuming uniform temperature. The coolant is assumed 
to have a constant, specified temperature, and the skin is 
heated convectively. A strong time variation of the con- 
vection coefficient, (Fig. 7b) simulates the passage of a 
moving shock. The resulting transient temperature and 
bar compressive stress exhibit features to be expected in 
a more complex model. 

The transient temperature history (Fig. 8a) shows a 
rapid rise following the sudden increase in the convec- 
tive coefficient, and then a smooth decay after the large 
convective heating is removed. The convective cooling 
causes the temperature to return quickly to equilibrium. 
The initial rapid increase in temperature causes the bar 
to yield in compression very early in the response (Fig. 
8b). After the temperature begins to decay at t = 0.5s, 
the stress rapidly changes to tension hence as the tem- 
perature returns to equilibrium a large residual tensile 
stress remains in the bar. The tensile stress is induced 
because the initial high temperature causes compressive 
yielding. This yielding tends to permanently shorten the 
bar. However, the bar’s fixed end boundary conditions 
prohibit the bar from decreasing in length as the tem- 
perature decreases so that a tensile stress develops as the 
material cools. 

During the thermal-structural analyses an uncer- 
tainly arose concerning applicability of the hardness evo- 
lution equation, Eq. (18), at high temperature rates 
(about 5000 °C/sec) occurring in these analyses. At such 
high temperature rates this equation may predict some- 
what anomalous behavior of the material. This is illus- 
trated by the following example. Consider a specimen 
of superalloy D 1900 + 7/ / with properties presented in 
the appendix. If the virgin, load -free specimen is heated 
from 760°C to L060°C in 0.1 second, the temperature- 
dependent variable Z 2 drops from 2700 MPa to 1200 
MPa. However, the actual hardness variable Z 1 , which 
follows Zj only through the relatively slow thermal re- 
covery term (second term in Eq. (18)) reduces in this 
time only to 2696 MPa. This means that although the 
specimen has reached an elevated temperature, it has not 
softened yet — a somewhat surprising result. 

The need for further research in thermal history ef- 
fects on D 1900 + // / was noted in ref. 6. 

In the initial viscoplastic analysis a fixed time step 
of 0.001 s was used, and 1200 time steps were required 
for the analysis. The variable time step algorithm shown 
in Fig. 6 was then implemented, and the problem was 
resolved. Fig. 8c shows the history of the variable lime 
step and indicates that the new analysis required only 
213 steps — a substantial savings. The figure shows that 
in the “fiat” part of the stress response a large lime step 
was used, but near t = 0 and again at l = 0.5s when the 


stress is changing rapidly, small time steps are needed 
to capture the response accurately. 'This example shows 
that the adaptive time step algorithm is an important 
step towards producing a reliable, accurate solution. 

Convectively Cooled Structure 

A more realistic 2D model of a convectively cooled 
hypersonic structure is shown in Fig. 9. The model 
represents a segment of a convectively cooled structure 
such as a wall of a scramjct engine fuel injection strut. 
The finite element thermal model (Fig. 9a) includes: (1) 
conduction heat transfer in the aerodynamic skin, heat 
exchanger fins and and primary structure, (2) convec- 
tive heat transfer between the walls of the coolant pas- 
sage and coolant, (3) mass transport convection in the 
coolant which has an unknown bulk temperature, and 
(<l) surface radiation on the aerodynamic skin. The aero- 
dynamic skin is uniformly convectively heated over its 
length, and superimposed on the uniform heating is a 
local, intense heating simulating a transient shock. The 
transient heating is induced by the time-dependent con- 
vection coefficient shown in Fig. 9b. 

In the plane strain structural finite element model, 
the primary structure and aerodynamic skin have unit 
thickness, but the heat exchanger fins in the coolant pas- 
sage are approximately represented by a single fin with 
the total thickness of 0.060 in. The wall segment has 
fixed displacement boundary conditions at the left and 
right ends. The aerodynamic skin, exchanger fins and 
primary structure are also loaded by internal pressure 
(7.0 MPa = 1000 psi) in the coolant passage which is rel- 
atively large compared to external aerodynamic pressure 
(0.35 MPa = 50 psi). 

The model was analyzed first for steady aerodynamic 
heating and internal pressure operating conditions. Then 
the sudden localized aerodynamic heating was applied. 
Thus, the steady temperatures and stresses serve as ini- 
tial conditions for the transient thcrmo-viscoplastic anal- 
ysis. 

The viscoplastic solution was computed with the vari- 
able time step algorithm, and the time step was varied 
using Eq. (39). A total of 265 steps were required to 
compute the response for a total time duration of 1.2 s. 

The thermal response of the wall segment is shown in 
Fig. 10. Figure 10a shows the lime history of the temper- 
ature at a point on the aerodynamic skin directly under 
the transient heating; Fig. 10b shows contours of temper- 
atures at / = 0.5s when temperatures are maximum. The 
temperature history is qualitatively similar to the results 
obtained in the ID model showing the rapid skin temper- 
ature rise and fall with the simulated shock heating. The 
temperature contours show relatively steep thermal gra- 
dients induced by the local healing, and that the coolant 
substantially limits the extent of the induced high tem- 
peratures. Thus the high temperatures are confined to 
the aerodynamic skin, and the primary structure experi- 
ences only small temperature changes. The temperature 
gradients in the skin particularly at the coolant-skin in- 
terface are not predicted with high accuracy because of 
the engineering model of the coolant heat transfer. How- 
ever, local temperature levels in the skin are reasonably 
accurate since net energy transfer to the coolant is mod- 
eled satisfactorily. 


Histories of the horizontal stress component a x at two 
points through the skin thickness arc shown in Fig. II. 
The stress histories follow the temperature, and under 
the intense local heating stresses are very similar to the 
results obtained from the one-dimensional model. At this 
location the skin yields through most of its thickness. Af- 
ter the heating ceases there is a rapid decay of stress, and 
under the intense local heating there are residual tensile 
stresses. Figure 12 shows the lime history of the vertical 
component of stress a y in the heat exchanger fins. High 
tensile stresses are induced with significant local yielding 
which are followed by residual compressive stress. The 
crude finite element model of the fins only approximates 
these stresses, but the high tensile stresses can poten- 
tially cause a bond failure at the heat exchanger/skin 
joint. The consequences of such a failure are studied in 
the next example. 

The viscoplastic stress histories are compared with 
stresses predicted assuming elastic behavior in Figs. 
13-14. The elastic computations were made with 
temperature-dependent elastic properties. The loss of 
stillness due to the elevated temperatures for 0 < I < 0.5s 
accounts for the stress “rolling-over” with an appear- 
ance of yielding. Generally, the elastic stresses are too 
high and, of course, since yielding does not occur residual 
stresses are not predicted. Maximum deformations (not 
shown) predicted by the clastic analysis and viscoplastic 
analysis are about the same, approximately 0.001 in. 

The two-dimensional character of the stress compo- 
nents a x and <t, are shown in Figs. 15 and 16, respec- 
tively. The stress distributions are shown at three times 
l — 0.05s, 0.5s, 1.2s in the response. The stresses at 
t = 1.2s are the residual stresses. Figure 17 shows con- 
tours of the principal plastic strains at t = 0.5s. The 
contours show the relatively localized nature of the high 
stress gradients as well as the tensile and compressive 
regions. The significant residual stresses suggest the pos- 
sibility of cumulative damage under repeated load cycles. 
This possibility needs further investigation. 

Convectively Cooled Structure With Damage 

The preceding analysis showed high tensile stresses 
in the heat exchanger fins and suggested the possibility 
of an aerodynamic skin/heat exchanger bond failure. To 
simulate such a possibility, the analysis was repeated but 
with the fins given greatly reduced material properties 
for a width of 0.30 in. under the high localized heating. 
Two results of this analysis are shown in Figs. 18-19. 

Since the fins can no longer support tensile stress, 
the internal coolant pressure causes substantial local de- 
formation of the aerodynamic skin. Figure 18 shows the 
deformed structure at t = 0.5s. The deformations, shown 
uuinagnified, represent a significant permanent displace- 
ment. This result can be seen in Fig. 19 which compares 
displacement histories for the undamaged and damaged 
structure. For the damaged heat exchanger case, a resid- 
ual permanent deformation of 0.007 in. is introduced by 
the thermal loading. 

The maximum deformation induced by the fin dam- 
age depends strongly on the width of the damaged re- 
gion. Increasing the width of the damaged region by a 
small amount will substantially increase the deformation. 


'['lie analysis of those large deformations requires that the 
present formulation be modified to account for large de- 
formation eirects. 

A possible consequence of this permanent deforma- 
tion is to cause a disturbance of the external flow. Pro- 
tuberances into high speed flows can cause significant lo- 
cal augmentation of the aerodynamic heating. Thus, the 
significant permanent deformation introduced by a fin 
bond failure is likely to cause a complex flow-tliermal- 
slructural interaction. 

CONCLUDING REMARKS 

A thcrmo-viscoplastic computational method for hy- 
personic structures subjected to severe local heating is 
presented. The method employs the Bodner-Partom 
unified viscoplastic constitutive model implemented in 
a finite element approach for quasi-static thermal - 
structural analysis. A variable time step algorithm is 
used to provide an efficient solution scheme for the stiff 
ordinary differential equation that characterize the evo- 
lution of the internal slate variables. 

An analysis of a simplified model of a one-dimension- 
al convectively cooled structure shows the basic features 
of the thermo-viscoplastic response and demonstrates 
the substantial computer time savings as well as the re- 
liability of the variable time step algorithm. An analy- 
sis of a two-dimensional model of a realistic convectively 
cooled structure provided more detailed understanding of 
the thermal structural behavior. The coolant flow dom- 
inates the thermal response providing a relatively short 
thermal transient. Under intense heating, significant lo- 
cal plasticity occurs in the aerodynamic skin, but the 
primary structure remains undamaged. Heat exchanger 
fins experience high tensile stresses and both the aerody- 
namic skill and heat exchanger fins have significant resid- 
ual stresses. An analysis of an aerodynamic skin/heat 
exchanger fin bond failure showed that the aerodynamic 
skin would experience significant local plastic deforma- 
tion due to the internal coolant passage pressure. The 
local deformation is pronounced enough to disturb the 
external aerodynamic flow and introduce a complex (low- 
llicrmal-structural interaction. 

The thermo-viscoplastic analysis of convectively 
cooled hypersonic structures under intense local heating 
shows the important role that such analyses can make 
in understanding complex transient inelastic structural 
behavior at elevated temperatures. 
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APPENDIX 

THERMAL-STRUCTURAL DATA 
Thermal Data 

1. Skin, heat exchanger fins and primary structure 
p = 0.283 lb m /in 3 

Temperature Dependent Conductivity and 
Specific Heat 


T 

K 

Cp 

l°«! 

[DTU/in - sR j 

[DTU/lb m - R\ 

0 

1.23 x 10" 4 

0.12 

500 

1.8 x lO” 4 

0.135 

1500 

3.25 x lO' 4 

0.177 

3000 

6.0 x 10" 4 

0.26 


2. Coolant 


m 

= 7.16 xlO' 3 

lb m / sec 

c P 

= 3.48 

DTU/lb m - R 

h 

= 665 

R 

h 

= 6.2 x 10- 3 

DTU/s - in J ■ 


3. Surface convection 

T x = 5345 R 

h = 6.0 x 10~ 5 DTU/s — in 7 - R 
Additional li at time 0 < t < 0.5 sec: 

0 < x < 0.375 : h aii = 0 

0.375 < x < 0.625 : h aii = 3.0 x I0~ 3 cos - ‘ 5) - 
~ " 0.05 2 

0.625 < x < 1.0 : h a u = 0 

4. Surface radiation 

<r = 3.30633 x 10' 5 BTU/s - in 2 - It* 
e = 0.8 
a = 0.8 


Structural Data 

(su 

peralloy U 1900 + 

1 . Temperature- 

indt 

-pendent constants 

Do 

= 

1.0 x 10 4 

sec -1 

z. 

= 

3000 

MPa 


= 

1150 

MPa 

>i 

= 

2 


rj 

= 

2 


m, 

= 

0.270 

MPa-' 


m 2 = 0 


2. Temperature-dependent parameters 


T 

\°C] 

< 760 

871 

982 

1093 

n 


1.055 

1.03 

0.850 

0.70 

II 

N 

[MPa\ 

2700 

2400 

1900 

1200 

d, 

(sec -1 ] 

0 

0.0055 

0.02 

0.25 

A* 

[sec -1 ] 

0 

0 

0 

0 





3. Elastic constants 


E = 1.987 x 10 5 + 16.78T - 0.1034T 1 + 
1.143 x 10' S T 3 [MPa] 

G = 8.650 x 10 4 - 17.58T + 0.02321T 3 
-3.464 x 10~ 5 T 3 [MPa] 

(with temperature in °C) 

4. Thermal expansion coefficient 


T 

a 

l'C\ 

1/*C]-10- 5 

0 

1.15 

200 

1.34 

400 

1.36 

600 

1.41 

800 

1.54 

1000 

1.60 

1200 

1.66 




Fig. 1: Fluid-thermal-structural interactions on an 
aerospace plane scramjet engine leading edge. 


Fig. 2: Convectively cooled structure. 
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(a) thermal-structural model 


(b) convection coefficient history 


Fig. 7: 1-D Model of convectively cooled structure 
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Fig. 8: Thermo- viscoplaslic response of I D model 
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(a) thermal-structural model 



Fig. 9: 2-D model of convectively cooled structure 
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Fig. 10: Thermal response of 2-D model 





Fig. 11: Viscoplastic stress histories in aerodynamic 
skin for 2-D model. 


Fig. 13: Comparison of viscoplastic and elastic stress 
histories in aerodynamic skin. 
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Fig. 12: Viscoplastic stress history in heat exchanger 
fins for 2-D model. 


Fig. 14: Comparison of viscoplastic and elastic stress 
histories in heat exchange^ fins. 
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Fig. 15: Contours of a x stress in 2-D model 
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(b) t = 0.5s 



Fig. 16: Contours of cr y stress in 2-D model 









Fig. 17: Principal plastic strain in 2-D model at t = 0.5s. 



l'ig. 18: Maximum dcfoi malion of slinclmo (t = 0.5s) with damaged heal exchanger fins. 
Deformation is shown to scale. 


u = 7.6*- 10‘ 5 in 

"** with fin damsge 



Fig. 19: Viscoplastic deformation of aerodynamic 
skin for undamaged and damaged heat exchange fins. 











